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1 INTRODUCTION 



Detailed studies on strong-interaction physics and the evaluation of the background for new-physics signals 
at high-energy colliders require accurate calculations in perturbative QCD. In this talk I shall describe 
some recent progress || in the study of a particular class of higher-order QCD contributions to hadronic 
collisions. This class concerns multiple emission of soft giuons. 

Some basic features of soft-gluon rcsummation are recalled in Sects. ^,|[ These Sections provide a 
brief introduction to the theoretical analysis of Ref. which is summarised in Sect. |[ Phcnomenological 
implications are reviewed in Sect, ^and some comments are left to Sect. ^ 

2 SOFT-GLUON EFFECTS 

The finite energy resolution of any particle detector implies that physical cross sections are always inclusive 
over arbitrarily-soft particles produced in the final state. This inclusiveness is essential in QCD calculations. 
Higher-order perturbative contributions due to virtual giuons are infrared divergent: the divergences are 
exactly cancelled by radiation of undetected real giuons. 

More precisely, if 1 — z denotes the fraction of the centre-of-mass energy y/S carried by unobserved 
final-state particles, virtual (v) and real (r) soft giuons affect the cross section by the following emission 
probabilities 

dwy{z) f^~' dz' 1 

- -2a(5(l-z) / -In- -, 



dz Jq 1 — z 

dwrjz) 1 1 

__ = 2a^ln^e(l-z-.), (1) 

where a — Cas/TT and the coefhcient C depends on the process. These double- logarithmic (DL) expressions 
arise from the combination of the customary bremsstrahlung spectrum duj/uj with the spectrum dO'^ for 
collinear radiation. Here I have introduced an unphysical cutoff e on the minimum energy fraction of both 
giuons. Adding the real and virtual terms, the physical limit e — )■ can be safely performed, leading to a 
finite differential probability: 



dw{z) 

— ; — = Za 

dz 



(2) 



where, as usual, the notation [(7(z)]-|_ stands for dz f{z) [g{z)\j^ = dz[f{z) — f{l)]g(z) , which defines 
a well-behaved distribution acting on any smooth function f{z) (such as cross sections or parton densities) 
at z = 1. 

Note that the virtual term contributes only at z — I while the real one, besides regularizing the virtual 
probability at z = 1, is spread out up to the kinematical boundary as given by the energy fraction x of 
the tagged final state. Thus the total soft-gluon contribution to the cross section is proportional to the 
following quantity 

'd^Mf)=_aln2(l_a:) . (3) 
dz 

This perturbative correction is the finite heritage of the cancellation of the infrared singularity. 

If we force the tagged final state to carry most of the total energy (i.e. we consider the quasi-elastic limit 
a; — > 1), we strongly suppress the radiative tail of the real emission. Thus the virtual term is unbalanced 
and the double- logarithmic factor in Eq. (||) can become large, asln^(l — a;)<;l, even if the coupling 
constant is in the perturbative region 1. In this kinematic regime, multiple soft-gluon radiation 

produces higher-order contributions of the type 

C„™ agln"(l-:E) , m < 2n , (4) 

that have to be resummed in order to obtain reliable theoretical predictions. 

In the case of hadron collisions, this large- a; region is encountered in the production of systems of 
high mass near threshold. Outstanding examples of these systems are the hadronic final state in 
deep-inelastic lepton- hadron scattering (DIS) (here x — xb is the Bjorken variable), lepton pairs with large 
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invariant mass produced via the Drell-Yan (DY) mechanism {x = t = Q'^/S) heavy quark-antiquark 
pairs (x = p — Attiq/S) Q and pairs of jets at large transverse momentum {x = t — Mj^/S) [Q. 

Similar soft-gluon effects occur both in hadron and lepton collisions when one measures detailed prop- 
erties ((5_L-distributions (|], event shapes 1^, multijet rates 0) of the hadronic final state. 

Although the logarithmically-enhanced contributions in Eqs. (HJ) are expected to be very relevant 
when x — !■ 1, the actual size of the soft-gluon corrections in cross section calculations depends on the 
coefficients Cnm and on the x-shape of the parton densities. Thus soft-gluon effects can be substantial also 
before reaching this kinematic region. This is the main motivation in Rcfs. |j for evaluating soft-gluon 
effects for top quark production at Tevatron |l0[ |ll|, ^ where the inelasticity variable p — Arri^/S is as 
small as jO = 0.04. 



3 RESUMMATION AND EXPONENTIATION 

In order to discuss soft-gluon resummation in hadron collisions, let me consider the DY process, which is 
the most studied process (see |^ and references therein) so far. The DY cross section <t{Q^, t) is obtained 
by convoluting the parton densities /(a;.;, Q^) of the colliding hadrons with a partonic cross section o'^Q^ , z) 
as follows 

1 

2a 



a{Q',T)= dz dxidx2fixi,Q')fix^,Q')Six,X2Z^r)A{z,asiQ'))MQ") ■ (5) 
Jo Jo 

Here, I have omitted all parton indices and factorized the partonic cross section in the Born-level contribu- 
tion (To times the term A that takes into account all the radiative corrections. The latter is computable in 
QCD perturbation theory as a series in as- In Eq. I have also explicitly introduced energy conservation: 
the lepton pair carries the fraction z = Q"^ / {X1X2S) of the centre-of-mass energy squared X1X2S available 
in the partonic subprocess. 

If the lepton pair is produced close to the hadronic threshold (t — s- 1), energy conservation forces 
the partonic subprocess towards z = \ (as well as the parton densities towards Xi = 1) and the radia- 
tive corrections in A(z,q;s) are dominated at lowest order by the soft-gluon probability in Eq. (^. The 
resummation of higher-order soft-gluon effects is conveniently carried out in A^-moment space. One can 
introduce the the A^-moments crN{Q^) of the cross section in Eq. (^ (and likewise for any other function 
of longitudinal- momentum fractions) by performing a Mellin transformation at fixed Q^: 

aN{Q')= f' drr^-' a{Q',T) . (6) 

JQ 

Working in A^-moment space, any z-distribution is replaced by a function of A^ and, due to the weight 
factor r , the threshold region is sampled by the limit A^ —> 00. However, the main reason for using 

A^-moments is not the replacement of distributions by functions, but rather the fact that one can easily 

implement momentum conservation. For instance, the energy-conservation constraint in Eq. is exactly 

diagonalized and the A^-moments of the cross section have the following factorized expression 

^jv(g') = fN{Q^)fN{Q^)AN{as{Q^)) MQ^) ■ (7) 

This kinematical factorization plays an essential role in the resummation programme. The radiative 
factor A{z,as) is obtained as follows 

+00 „i 



A(z, as) = (5(1 - z) + / dzi. .. dzn '""^ ^' j ' 6ps(2; zi,.. z„) , (8) 
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where the probability dwn for producing n soft gluons is integrated over the phase-space region (symbolically 
denoted by 8ps) that is available according to the actual definition of the cross section. The first relevant 
ingredient to perform the summation in Eq. (j8|) has a dynamical origin. Owing to the factorization 
properties of multi-giuon QCD amplitudes in the soft limit, the n-gluon probability factorizes in the product 
of the single-gluon contributions dw{zi) in Eq. (||): 

dWnizi, . . . , Zn) I -fj dw{zi) 



dzi . . . dzn n! 
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The second relevant information for the summation regards the kinematics. In general, the phase-space 
function Ops depends in a non-trivial way on the multi-gluon configuration so that it cannot be handled 
in a simple manner. However, in the case of total cross sections, like in the DY process, the only relevant 
constraint is longitudinal-momentum conservation. This constraint is exactly factorizable going to A^-space: 

Qps{z; Zi, . . . , Zn) ^ S{z - Zi . . . Zn) — > z^^^...z^^^ (10) 

A^— moments 

Using Eqs. (^ and ( px| ) one can straightforwardly recast Eq. (^) in exponential form: 



A AT (as) = exp 



dz z 







dz 



(11) 



In summary, the physical basis for the exponentiation of soft-gluon contributions are the factorization of 
both the multi-gluon amplitudes and the phase-space. The first property follows from QCD dynamics and 
is completely general. The second property depends on the actual definition of the cross section and can 
be easily violated even in QED jl^. In the case of total cross sections in hadron collisions, because of 
energy conservation, phase-space factorization can only be achieved by working in A^-space. There are no 
physical basis for exponentiation directly in x-space. 

Inserting Eq. (|^) into Eq. ( pi] ) and performing the large- A'' limit, one finds the exponential of a DL 
expression: 



^Nias) = exp 



^ z^-i - 1 1 
2a I dz In ■ 







(12) 



1 - z 1 - z_ 

= exp [-a (ln2 7V + 0(lniV))] . (13) 

The derivation of Eq. (|l^) that I have sketched would be complete in massless QED (neglecting the effect 
of soft-fermion pairs). The only simplification I have introduced in the QCD case regards the factorization 
formula (^). Unlike photons, which have no electric charge, giuons carry colour charge and multiple 
gluon emission is affected by dynamical correlations. Nonetheless, using general properties related to 
gauge invariance and unitarity, one can prove that, in the DY process and similar total cross sections, a 
generalized exponentiation theorem is valid in the following form 

lnAA,(as) = In A^. 91 (as In iV)+ ^2 (as In TV) + as 53 (as In TV) -f... . (14) 

Since the functions gi depend only on as In A^, in this context exponentiation means that all the logarithmic 
corrections of the type ag In™ A^ with n + 1 < m < 2n, allowed by Eq. (^, are simply taken into account 
by the exponentiation of lowest-order terms and thus cancel in the exponent of A^v- Owing to the general 
structure of Eq. ( p^ , one can control the resummation of the logarithmically-enhanced terms by introducing 
an improved perturbative expansion. The function gi gives the leading logarithmic (LL) contributions 
ag In""*"^ A^, g2 contains the next-to-leading logarithmic (NLL) terms a^ In" A^, etc. 

The functions gi essentially represent the QCD correction to the factorization in Eq. (||). Most of the 
correction is simply taken into account by introducing running coupling effects. Roughly speaking, this 
amounts to the replacement a = Cas((5^)/7r — > Cas((l — z)Q^)/7r at the integrand level in Eq. (p^). 
Thus, because of the logarithmic behaviour as(/i^) ^ 1/ ln(/i^/A^), integrating exactly in z up to z = 1 
one eventually hits the Landau pole at z = 1 — This singularity, called infrared renormalon [p^, p^, 

signals the onset of non-perturbative phenomena and has to be properly regularized. 



4 IMPLEMENTING RESUMMED FORMULAE IN HADRON 
COLLISIONS: THE MINIMAL PRESCRIPTION 

The theoretical analysis in Ref. deals with difficulties that arise when one tries to apply resummed 
formulae of the type discussed in Sect. ^ to the actual evaluation of physical cross sections. 

The first difficulty is due to the fact that these formulae are provided in A^-space. Cross sections 
like that in Eq. (^) are not experimentally measurable in the entire kinematic range < r < 1 and thus 
theoretical predictions for the A^-moments cannot directly be compared with data. In any phenomenological 
application one has to go back to x-space and compute the radiative factor A(z, as)- 
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The regularization of the running couphng singularity noticed at the end of Sec. ^ leads to the second 
difficulty. One has to introduce a regularization prescription that does not spoil the main perturbative 
features. 

As discussed in Ref . , several methods used in the literature to solve these difficulties are theoretically 
unjustified and, typically, enhance soft-gluon effects long before the hadronic-threshold region is actually 
approached. 

In particular, we pointed out the differences between the above difficulties and we warned of the danger 
of exploiting the formal equivalence N ~ 1/(1 — x) to straightforwardly transform A^-space resummed 
formulae into x-space resummed formulae. For instance, following Refs. @, ^, D, one would replace the DL 
expression (|l^) with the x-space version: 

A(z, as) ~ [e(l -z-e) e-'^'"'^!--' (1 + • ■ •)] , (15) 

where the dots stand for 'subleading' logarithmic terms, i.e. terms at most of the order of agln"(l 



z 



Expanding Eq. (15) as a power series in as and identifying e with the unphysical cutoff in Eq. (|l|), 
one obtains the lowest-order probability in Eq. (^. Equation (|l^) is formally valid ||] in terms of the 
logarithmic expansion in ln(l — z) , but the identification In In 1/(1 — z) and the ensuing approximation 
of neglecting subleading logs is correct only if In 1/(1 — z) ^ 1. Care has to be taken before extending 
this approximation to higher values of z. 

In particular, it is quite dangerous to insert Eq. ( p^ into the factorization formula (^, where the z- 
integration range extends up to z = 1. In the case of hadron collisions, after having introduced the parton 
densities f{xi,Q^) in customary factorization schemes (like the DIS or MS schemes), the coefficient C in 
a = Cas/TT turns out to be negative. Thus the expression exp[— aln^(l — z)] = exp[|a| ln^(l — z)] diverges 
faster than any power of (1 — z) and is not integrable at z = 1: in the resummed expression ( p^ the 
unphysical cutoff e can no longer be removed without producing a singularity. Note that such a singularity 
is not related to the hadronic threshold. The inelasticity variable r controls the lower limit of the z-integral 
in Eq. (^): the point z = 1 is inside the integration region even asymptotically far (r — )■ 0) from threshold. 

The origin of this singularity was studied in detail in Ref. Q] (see also [^) and traced back to the 
same-sign factorial divergence of the 'subleading' terms systematically neglected in the derivation of the 
x-space formula (^5|) from Eq. (^3|). No analogous divergence is present in the original A^-space formula. 
Soft-gluon exponentiation in A^-space is related to momentum conservation. Using exponentiation directly 
in x-space, one violates infinitely many times momentum conservation and can build up a (kinematical) 
divergence. 

In Ref. Q we specified a prescription for the implementation of soft-gluon resummation formulae in the 
computation of cross sections in hadronic collisions. We used the generalized exponentiation formula in 
Eq. (|l^) and performed numerically the inversion from A^-space to cc-space without further approximations. 
As for the regularization of the Landau pole in the running coupling, we pointed out that the functions 
gi{oLs In A^) are well-defined analytic functions of N in the complex plane. The Landau pole only produces 
a branch point at A^ = Q/A. In the numerical inversion we used a prescription that leaves this non- 
perturbative branch cut to the right of the integration contour. 

As proved in Ref. Q, with our prescription no factorial divergence is introduced in the resummed 
expansion. Factorially growing terms related to infrared renormalons are likely to be present in the full 
perturbative QCD series, but as shown in Ref. ]T^ , their identification and evaluation require an analysis 
whose accuracy goes beyond that used to prove the exponentiation of logarithmically-enhanced contri- 
butions. Since our prescription fulfils the relevant kinematical constraints and does not introduce large 
corrections that are not justified by the logarithmic expansion, we call it the 'Minimal Prescription' (MP). 
This prescription can be used to extend perturbative QCD calculations towards the threshold region by 
consistently taking into account logarithmically-enhanced soft-gluon effects. 



5 PHENOMENOLOGICAL STUDIES 

In Ref. [Q, the MP was applied to the DY process and to the production of heavy quarks and high-p^ 
jets in hadron collisions. Since the last two processes are becoming increasingly topical, in this Section I 
briefly review our results. In both cases, in the resummation of soft-gluon corrections we included only 
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Figure 1: Contribution of gluon rcsummation at order and higher, relative to the NLO result, for the 
individual subprocesses and for the total, as a function of the top mass in collisions at 1.8 TeV. 



the LL contributions (i.e. the analogue of the function in Eq. (p^). In the case of heavy-quarks, this 
corresponds (roughly speaking) to the iV-space version of the resummation considered in Ref. Q| . 

5.1 Heavy-quark production 

The importance of the resummation effects for top quark production in pp-coUisions was studied by com- 
puting the following quantities 

'^gg '^qq '^gg + '^qq j-j^g-j 



Aaa) ' (qq) ' Agg) , Jqq) 

'NLO "NLO "NLO ~^ "NLO 



that are plotted in Fig. 1. Here Sat is equal to our MP-resummed hadronic cross section in which the terms 
and q;| have been subtracted, and 



of order q;| and q;| have been subtracted, and ctI'^tci-, is the hadronic cross section in full next-to- leading 



order ||17| (i.e. including both soft and hard radiation up to order Og). The labels ab refer to the various 
partonic channel that contribute to the cross sections. 

At Tevatron energies one can see that the contribution of resummation is very small, being of the order 
of 1% for mt — 175GeV (in the calculation of Ref. jo) the resummation effect is about 10%, i.e. one order of 
magnitude higher; similar quantitative effects are found with the method of Refs. D). Thus, as expected 
by the small value of p = Arri^/S ~ 0.04, we concluded that top quark production at Tevatron is reliably 
estimated by the next-to-leading order (NLO) QCD calculation. Based upon these findings, we updated 
the computation of Ref. (isj . Full details of our calculation are given in Ref. |^ . The QCD prediction is 
conveniently parametrized as follows 



175- 



atf(1.8TeV) = e^r^(4.75:^;^^)pb. (17) 

We evaluated others heavy-quark production cross sections. In general we found § that, in most 
experimental configurations of practical interest, soft-gluon resummation effects are not dominant either 
because they are very small or because they are well below the (estimated) uncertainty due to higher-order 
(non-soft) corrections. 



5.2 Jet production 

The interest in the effects of soft-gluon resummation on the behaviour of jet cross sections is prompted 
by the discrepancy between the single-inclusive jet distribution at large Pt, as measured by CDF and 
the result of the NLO QCD predictions This discrepancy is the topic of ongoing investigations and 
discussions ||2^ on both experimental and theoretical aspects. 

Owing to its kinematic similarity with DY and heavy-flavour production, we studied [Q the soft-gluon 
corrections to the invariant-mass distribution of a jet pair. For this distribution an analogous discrepancy 
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between data and theory has been reported |^ . Figure 2 shows our results for the following quantities 

.(4) c(4) c(4) X(4) , r(4) .(4) 

Qgg "qg "qq Ogg ^ "qg ^ "qq , . 

ct(3) ' ct(3) ' o-(3) ' ct(3) ' ^ 

where, analogously to Eq. (^6|), S^^'' is equal to the MP resummed hadronic cross section with terms of 
order a| subtracted. Unlike in Eq. (|l6|), ct*^'^^ is not the full NLO cross section but rather its approximation 
as obtained by truncating the resummation formula at order a^. One can see that the resummation effect 
leads to 10 15% increase of the cross section when the inelasticity variable r — Mjj/S is as large as 
T - 0.5. 

These results should only be taken as an indication of the order of magnitude of the correction. This is 
because we did not include a study of the resummation effects on the determination of the parton densities, 
we only considered LL resummation and we did not implement detailed experimental cuts. 

Note also that other distributions, such as the Pt of the jet, have a rather different structure from 
the viewpoint of soft-gluon resummation because of additional jet-broadening effects due to final-state 
radiation. Calculations of these effects are in progress. 

6 SUMMARY AND OUTLOOK 

The MP can be used to evaluate the effect of soft-gluon resummation on cross sections in hadron collisions. 
Using the MP we find that the logarithmically-enhanced soft-gluon terms become important only fairly 
close to the hadronic threshold (x ~ 0.5). Further investigations of resummation effects in this kinematic 
region are warranted. These include the consistent determination and Q^-evolution of the parton densities 
by using resummed anomalous dimensions Q and the evaluation of NLL corrections ||2^ . 
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